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Abstract 



An ordinary differential equation is said to have a superposition 
formula if its general solution can be expressed as a function of a finite 
number of particular solution. Nonlinear ODE's with superposition 
formulas include matrix Riccati equations. Here we shall describe 
discretizations of Riccati equations that preserve the superposition 
formulas. The approach is general enough to include g-derivatives 
and standard discrete derivatives. 
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1 Introduction 

The well-known Riccati equation 



arises in numerous contexts. For example, one can find particular 
cases of (fQ) as equations describing Darboux transformation for the 
Schrodinger operator X, as equations describing particular solutions 
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w = a(t) + b(t)w + c(t)w 2 




of Painleve equations [2], or as Backlund transformations for soliton 
like equations [31 E] ■ 

The Riccati equation has many interesting properties, one of them 
is the fact that it is a non-linear equation with a superposition formula. 
A differential equation has a superposition formula if it is possible 
to express its general solution as a function of a finite number of 
particular solutions and arbitrary constants. 

It is easy to prove that the Riccati equation Q has a superposi- 
tion formula. Let wo,wi,W2 be three different particular solutions of 
Eq. (fT]l. It is well-known that using one particular solution one can 
transform Q into a linear equation. Indeed, let w = y + Wq, then we 
have y = (b + 2cwo)y + cy 2 . Now let y = j, then we obtain a linear 
equation x = —(b + 2cwq)x — c. We know two particular solutions 
%i = 77, 1 „ , and xo = ; , -; , of this linear equation, hence we can 
find a general solution of this linear equation using xi, xi and one 
arbitrary constant. This gives us a formula for the general solution of 
the Riccati equation in terms of three particular solutions wq,w\,W2 
and one arbitrary constant c : 

{w 2 - w ){wi - w ) 

W = W + —, r. (2) 

W2 - W + C{Wl - W2) 

One can also consider the matrix Riccati equation 

W = A(t) + B(t)W + WC(t) + WD(t)W, (3) 

where W is a n x /c-matrix and ^4, B, C, D are matrices of appropriate 
sizes. This equation also has a superposition formula [31 E]. 

One can discretize the Riccati equations (^Q) and (J5J) in different 
ways, but we are interested in discretizations possessing superposition 
formulas. Our main result is the following. Let us consider a two- 
parameter class of operators U q ^ : 

[ f'{t) if q = 1, h = 

U q ,hf(t) = < ffqt +h )-f(t) . (4) 

y ' 1 (^zTjTt+TT 111 other cases. 

As we shall show below the equation 

U q , h w{t) = a(t) + b(t)w(t) + w(gf + /i)c(t) + w(gi + h)d(t)w(t) (5) 

has a superposition formula. In the case where q = l,h = we 
obtain the Riccati equation (^Q). It can be easily seen that in the case 
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q = 1, h ^ we obtain the /i-Riccati equation 



w(t + h) 



w(t) 



a(t) + b{t)w(t) + w(i + h)c{t) + + h)d(t)w(t) 



h 



and in the case h = 0,q ^ 1 we obtain the g-Riccati equation 



Analogous results hold also in the matrix case. 

It is interesting to remark that particular cases of such discretiza- 
tions appear in the same contexts as their classical analogue, for ex- 
ample as Darboux transformations for the discrete Schrodinger oper- 
ator [D , or as equations describing particular solutions of discrete 
Painleve equations [SJ. Thus the discretizations constructed in the 
present paper are very natural. 



2 Differential and difference equations 
with superposition formulas 



Let us consider a (vector) first order ordinary differential equation 



We shall say that this equation has a superposition formula if its 
general solution y(t) can be expressed as a function of a finite number 
m of particular solutions y i , . . . , y m and n free constants 



The formula ((JJ) is called a superposition formula for Eq. (|5|>. 

The study of ordinary differential equations with superposition 
formulas has a long history that goes back to Sophus Lie. He proved 
the following fundamental theorem. 

Theorem 1 (Lie |9j) The equation Jj^ has a superposition formula 
if and only if the function f has the form 



w(qt) — w(t) 
(9-1)* 



a(t) + b(t)w(t) + w(qt)c{t) + w(qt)d(t)w(t). 



y k = f k (y,t), k = l,... 



n. 



(6) 




(7) 



r 



/ fc (y,i) = J>(i)£f(y) 



i=i 
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where the functions £z(y) are such that the vector fields 

n 
k=l 

generate a finite dimensional Lie subalgebra f) of the algebra of vector 
fields ofC N orR N , i.e. 

r 

[Xi,Xj] = ^2 CijkXfr. 
fc=l 

Difference equations with superposition formulas can be denned in 
a similar way. In order to see what can be done in the difference case, 
let us reformulate the Lie theorem in different terms. 

Let G be a Lie group acting on a manifold M. Let us consider a 
curve g(t) on G such that g (to) = e and a point Uq € M. It is easy to 
find a differential equation for a curve u(t) = g(t) • Uq on M. Indeed, 

u(t) = g(t) • n = <?(%W~Vi) • u = g{t)g(t)- x ■ u(t). 

Since gg~ l £ 0, we see that the following statement holds. 
Statement 1. A curve u(t) = g{t) ■ uq, such that g(to) = e, is a 
solution of a differential equation 

ii(t) = £(t) ■ u(t), (8) 

where £ = g r 5'~ 1 S 0. The point uq plays the role of an initial condition, 

u{t ) = Uq. 

Let us now consider Eq. (jHJ) with given Let be a solution. 
Let us choose a fixed value of the parameter to and the initial value of 
the solution uq = u(to). There exists a unique curve g(t) S G such that 
gg~ x = £ and g(to) = e. We see that the curve u(t) = g(t) ■ uq is also 
a solution of Eq. (JHJ, and it also satisfies the initial condition u(to) = 
uq. The following statement follows from the uniqueness theorem for 
solutions of ODEs. 

Statement 2. Let £(t) be a curve on q. Let to De a fixed value 
of the parameter. Any solution u(t) of equation has the form 
u(t) = g(t) ■ no, where no € M is the initial condition u(to) = uq 
and g(t) is a curve on G such that g(to) = e. The curve g(t) S G is 
uniquely determined by the conditions gg~ l = £, and g(to) = e, hence 
g(t) is the same for all solutions u(t) of equation (jHJ). 
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Statements 1 and 2 mean that Eq. (JSJ) has a superposition for- 
mula. Indeed, let ui(t),... ,u m (t) be particular solutions of ©. Let 
us choose an initial value of the parameter to- It follows from State- 
ment 2 that ui(t), . . . , u m (t) have the form Ui(t) = g(t) ■ Ui{to), where 
g(t) is the same for all i. Then we have the following system of equa- 
tions for g(t) : 

u m {t) = g(t) ■ u m (to). 

If m is sufficiently large, g(t) can be expressed in terms of Ui{t) and 
Ui(t ) : 

g(t) = F(ui(t), . . . ,u m (t);ui(t ), . . . ,u m (t )). 

It was shown by Lie that m should at least satisfy the inequality 
mn > r, where n = dimM and r is the same r as in Theorem ^ i. e. 
the dimension of the image of q in the algebra of vector fields on M. 
It follows that the general solution has the form 

u(t) = g(t) ■ u = F(ui(t), . . . ,u m (t)]Ui(t ), . . .,u m (t )) ■ u(t ), (9) 

the initial condition u(to) plays the role of the arbitrary constant in 
the superposition formula (|9*|). It should be remarked that the function 
F in the superposition formula © does not depend on the choice of to, 
hence the superposition formula is essentially the same for any choice 
of to- The superposition formula obtained by choosing another initial 
value of the parameter t\ can be obtained form (JHJ) by expressing the 
constants U\(to), ■ ■ ■ , u m (to),u(to) in terms of iti(ti), . . . , u m (t\), u(ti). 

It should also be mentioned that generally we also have some inde- 
pendence conditions on the particular solutions in order to find g(t). 
For example in the case of a homogeneous linear equation of second 
order any solution is a linear combination of two particular solutions 
under the condition that these two solutions are linearly independent. 

Statements 1 and 2 are not complicated. However, the Lie theorem 
also contains the following statement which is more difficult to prove. 
Statement 3. Every ordinary differential equation with a superposi- 
tion formula must have the form (JHJ) for some group G acting on some 
manifold M jOJ. 

A modern exposition of the proof can be found in |lUj . State- 
ments 1, 2 and 3 are exactly what Lie's theorem states. There is a 
very important thing used implicitly. In the proof of Statement 2 
a crucial step makes use of the uniqueness theorem for solutions of 
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ODEs. We are working implicitly with some class of functions such 
that for solutions from this class this theorem holds. In the same im- 
plicit manner we assume that equation (jHJ) is an equation for which 
the uniqueness theorem holds. For differential equations the assump- 
tions underlying the uniqueness theorem are very natural. One can 
assume that Eq. (jSJ) is defined by a continuous curve and that we 
are considering solutions in the class C 1 (M), i. e. the class of differen- 
tiable functions with continuous derivative. For difference equations 
the situation is different and the appropriate class of solutions must 
be carefully chosen and we need a more rigorous definition. 
Definition. A differential (difference) equation is said to have a su- 
perposition formula in some class of solutions if the general solution 
belonging to this class can be expressed as a function of a finite num- 
ber m of particular solutions yi, . . . , y m , also belonging to this class, 
and free constants 



where n = dim M. Eq. (jlOj) is called a superposition formula. 

The question of describing all difference equations with superpo- 
sition principles is beyond the scope of this paper. It is a very dif- 
ficult question, especially since different difference operators require 
the consideration of completely different classes of solutions. We will 
see examples later in the case of the discrete Riccati equation. Let 
us concentrate on Statements 1 and 2 and see how one can construct 
difference equations with superposition formulas. 

We can start again from a Lie group G, a manifold M, a curve 
g(t) G M such that <?(io) = e and a point uq € M. Let u(t) = g(t) ■ uq. 
Let U be a first order difference operator, for example of the form Q 
defined in the Introduction. The difference derivative Uu(t) can be 
defined only in a fixed coordinate system on M. It is no longer an 
invariant object like the vector u(i). It is impossible, in general, to 
write something like Uu(t) = Ug(t) ■ uq and reduce all to the Lie 
algebra, since there is no general chain rule for a difference operator 
and Ug(t) is not, in general, an element of q. However, in some cases 
we can construct surrogates of statements 1 and 2 in the following 
way. Let us assume that for some G and M (with a fixed coordinate 
chart) the following statements are true. 

Statement 1'. A curve u{t) = g(t) ■ uq such that g(to) = e is a 
solution of some difference equation 



y = F(yi 




(10) 
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Statement 2'. Let to be a fixed value of the parameter. Any solution 
of equation belonging to some class of solutions, has the form 

u(t) = g(t) ■ u , where u £ M is the initial condition u(t ) = uq 
and g(t) is a curve on G such that g(to) = e. The curve g(t) € G is 
uniquely determined by the condition g(to) = e, i. e. g(t) is the same 
for all solutions u(t) of the equation 

In this case we can easily prove in the same way as for ODEs that 
the difference equation (jllj) has a superposition formula in the class of 
solutions mentioned in the Statement 2'. Moreover, this formula will 
be the same superposition formula ® as for the differential equation, 
since in both cases we are working with solutions of the same form 
u(t) = g(t) ■ u . 

Let us now show how this works in the case of the matrix Riccati 
equation. 

3 Differential and difference Riccati 
equations 

It is easier in this Section to consider the matrix Riccati equation from 
the beginning. We will omit the word "matrix" for simplification. 

Let us now recall how the Riccati equation arises from the con- 
struction described in Section [2 Consider G = GL{N) and let M be 
a Grassman manifold G n +k,k of ^-planes in an N = (n+A;)-dimensional 
space. The homogeneous coordinates of a point p € M are given by 
the components of an (n + k) x /c-dimensional matrix 



where X is a n x A;-matrix and Y is a k x /c-matrix. The columns 
of 1)12|) span a A;-plane defining p. The point p is thus identified with 



identifying different bases for the same /c-plane. The action of an 
element 




(12) 




under the relation 
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upon M is obtained by the projection 



7T 



X 
Y 



X 
Y 



from a linear action 
9 ■ 



X \ _ / M N \f X 
Y -\ P Q [y 



On the afhne subspace defined by detY ^ we may define inho- 
mogeneous coordinates W = XY^ 1 , W is an n x A;-matrix. In the 
inhomogeneous coordinates the group action is given by the formula 

g - W = (MW + N)(PW + Q)~ x . (13) 

Let us now consider Eq. (jHJ. In this case it reduces to Eq. (j3J), 

where A, B, C, D are defined by 



-D(t) -C(() 

Let us now construct a difference version of the Riccati equation. 
We use the two-parameter class of operators U q h defined by the for- 
mula (j3J) in the Introduction. This class of operators can be described 
in an axiomatic way by the following properties: 

1. U q>h l = Q. 

2. Ug jh (af + bg)(t) = aU qjh f(t) + bU qjh g(t), where a and b are 
constants. 

3. [U qth ,t]f(t)=f(qt+h). 

For us the important property of these operators is the modified Leib- 
nitz rule: 

U q , h (fg)(t) = g(t) U q , h f(t) + f(qt + h) U q , h g(t). 

Let us proceed as described in Section |21 Let G = GL(N) and M 
be the Grassman manifold G n+ k,k as before, where N = n+k. We need 
to choose a coordinate chart on M : let us take the inhomogeneous 
coordinates described above. 

Let g(t) be a curve on GL(N) such that g(to) = e and Wq be the 
inhomogeneous coordinates of a point on M. We will represent g(t) as 
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M(t) N(t) 
P(t) Q(t) 



Let us find a difference equation for W(t) = g(t) ■ Wq. The group 
action is described by the formula (|13f) . so we have 

w(t) = (M(t)Wb + N(t))(p(t)w + get))- 1 . 

We can find U q h W(t) directly, but it is better to proceed in the follow- 
ing way. Let us choose homogeneous coordinates for Wo, i. e. matrices 
X , Y such that XqY^ 1 = Wq. Then we can consider X(t) and Y(t) 
defined by 

( X(t) \ _ ( Xq\_( M(t)X + N(t)Y \ 
V Y(t) J 9{t) V Y Q J V P(t)X + Q(t)Y J ■ 

It is easy to see that W(t) = X(t)y(t) _1 . Since the group action 
in homogeneous coordinates is simply the matrix multiplication, it is 
easy to see that 

U q ,h ( y$ ) = (U q>h g(t)) ( ^ ) = {U q , h g{t))g{t)- 1 ( *® ) . 
Let us represent {U q ^g{t)) g{t)~ l in the following way: 

in the case of U\fi = -j- the matrix functions A(t), . . . ,D(t) are the 
same as in (|14|) . We obtain the following linear difference (or differ- 
ential in the case {7i,o) equations: 

U q>h X(t) = B(t)X(t)+A(t)Y(t), (15) 
U q , h Y(t) = -D(t)X(t)-C(t)Y(t). (16) 

It is easy to prove that for two matrices A(t) and B(t) we have the 
identity 

U q , h [AB- 1 ] (t) = [XJ q , h A{t)]B{t)- l -[AB- x ] {qt+h) [U q , h B(t)\ B(t)~ 

Using this identity and Eqs. (f 151 116|) we can finally find a difference 
equation for W(t) : 

U q , h W{t) = U q>h [XY- 1 ] (t) = 
= [U qth X(t)} Y (t)- 1 - [XY- 1 ] (qt + h) [U q>h Y{t)\ Yit)- 1 = 
= A(t) + B(t)W(t) + W(qt + h)C(t) + W(qt + h)D(t)W(t). 
Thus we have found Statement V from section [21 in our case: 
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Theorem 2 A function 

W(t) = (M(t)Wb + N(t))(P(t)W + QCt))- 1 , (17) 

such that 

and 5 (to) = I, is a solution of the difference matrix Riccati equation 

U q h W(t) = A{t) + B(t)W{t) + W{qt + h)C{t) + W(qt + &)!>(*) W(t) 

(19) 

with the initial condition W(to) = Wq. 

The matrices A(t), . . . ,D(t) are defined by the formula 



B(t) A(t) 
-D(t) -C(t) 



[U q , h g{t)]g{t)-\ (20) 



In the case of = -4- we obtain the usual (matrix) Riccati 
equation. 

It is easy to see that the limit of Eq. (|19j) when q — > 1 and h — > 
is the usual Riccati equation (j3J) with the same A(t), ... ,D(t). Hence 
we can consider Eq. (fl9|) as a discretization of the usual Riccati equa- 
tion ©. 

As explained in Section in order to prove the existence of a 
superposition formula we need to prove Statement 2' from SectionElin 
our case, i. e. we need to prove that any solution of Eq. lfTU|) belonging 
to some class of solutions has the form (|17j). where the point Wq is 
the initial condition W(to) = Wq, and the curve g(t) (fTH|) is uniquely 
determined by the condition g(to) = I, i. e. g(t) is the same for all 
solutions. 

As was already explained, we need a uniqueness theorem to prove 
this. This requires some assumptions on equation ()19|) and on the 
class of considered solutions. 



4 Superposition formulas 

The theory of superposition formulas for the discrete Riccati equa- 
tion (|19j) depends on the values of q and h. We will consider two very 
different cases: the case \q\ ^ 1 and the case q = l,h ^ 0. We will 
not consider other cases since the case q = 1, h = is simply that 
of differential equations, and the case \q\ = l,q ^ 1 would merit a 
separate treatment. 
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4.1 Superposition formulas in the case \q\ ^ 1 

The key observation in this case is the following. Let T be the trans- 
formation T '. t i — ► qt ~\~ h of the complex plane. Then the set of 
points T n (to),n € Z, has an accumulation point j^—. The well-known 
uniqueness theorem for holomorphic functions says that if a func- 
tion /, holomorphic in a disc, is equal to zero on a subset possessing 
an accumulation point inside the disc, then / is identically zero. This 
provides us with a uniqueness theorem for solutions of a difference 
equation in the class of functions meromorphic in the entire complex 
plane and holomorphic in a neighborhood of the accumulation point. 
Let us consider two solutions of Eq. (|19|) from this class. Let us con- 
sider a disc such that both solutions are holomorphic in this disc. For 
simplicity let us consider the case |g| < 1. Let our solutions be equal 
in some point to inside the disc. Then we shall prove that two solu- 
tions are equal also in the points T n (to),n G N. It follows from the 
above mentioned uniqueness theorem for holomorphic functions that 
two solutions are identical. 

On the other hand the class of functions meromorphic in the entire 
complex plane and holomorphic in a neighborhood of the accumulation 
point is very natural since if g(t) is holomorphic then the solution 
W(t) ()17|) is, in general, a meromorphic function. 

This leads us to the following theorem. 

Theorem 3 Let \q\ ^ 1. The discrete Riccati equation \19\) with holo- 
morphic coefficients A(t),B(t),C(t),D(t) has a superposition formula 
in the class of functions meromorphic in the entire complex plane and 
holomorphic in a neighborhood of the point 

Proof. Let us remark that the change of coordinates t = t' + j*— 
permits us to reduce Eq. (fTT?|) to an equation of the same form but 
with h = 0. Hence we shall consider only this case, the case of q- 
derivative. Let us also remark that the change of coordinates t = -t' 
permits us to reduce further to an equation of the same form but with 

M < i- 

Let us consider the formula (|20[) as a difference equation for the 
N x iV-matrix g(t). We can rewrite it as 

=W*)=(_ B 4) % )*>. (2D 

It is a homogeneous linear ^-difference equation for g(t). Let us also 
fix some to an d consider the initial problem g (to) = I- 
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It can be proven by methods similar to those used in Ref. ^2] that 
there exists a unique solution g(t) holomorphic in the entire complex 
plane. Let us prove that g(t) is non-degenerate, i. e. g(t) £ GL(N). 
Suppose that in some point t\ we have det g{t\) = 0. We can rewrite 
equation (|2*T]) as 



We see that if g(t) is degenerate then g(qt) is degenerate. It follows 
that for the holomorphic function det g(t) we have det g{q n ti) = 0, n € 
N. This implies g(t) = 0, but it contradicts the initial condition g(to) = 
I. Hence g(t) G GL{N). 

Let Wi(t) be a solution of Eq. (|19j) meromorphic in the entire 
complex plane and holomorphic in a neighborhood of the point 0. 

Let us choose some R > such that the following conditions hold 
in the disc \t\ < R : 

1. Wi(t) is holomorphic; 

2. The matrix 



is non-degenerate. 

Such R exists since C(t),D(t) and W\{t) are holomorphic in zero. 
Let us choose some to such that \to\ < R. Let Wq = W\(to) and let 



be the solution of (|21j) with the initial condition g(to) = I. Let us 
consider the function 



W 2 (t) = g(t) ■ W = (M(t)W + N(t))(P(t)W + Q(t))~\ (22) 



It follows from Theorem |2 that the function (|22[) satisfies Eq. (|19|) 
and the initial condition W^C^o) = Wo- Let us note that W%(t) is 
meromorphic in the entire complex plane. 
We can rewrite Eq. (|19j) in the form 



W(qt)[I-(q-l)t(C(t)+D(t)W(t))] = ( q -l) t (A(t)+B(t)W(t))+W(t). 
Since the matrix I — (q — l)t(C(t) + D(t)W(t)) is invertible in the 



Wi(gto) = [(q-l)t(A(to)+B(t )Wo)+Wo][I-(q-l)t(C{to)+D(to)W ^ 

12 




/-(?-l)t(C(t)+D(t)lfi(f)) 




disc \t\ < R and \q\ < 1, we have 



The same is true for W%(qt), so W2(qto) = W\{qto). We can prove in 
analogous way that W2(q n t) = Wi(q n t),n G N. 

Since W\{t) is holomorphic in a neighborhood of 0, we have 

lim W 2 {q n t ) = lim W^to) = Wi(0) 

n— >oo n— *oo 

and this implies that W 2 (t) cannot have a pole in 0. Since both W\(t) 
and W2(t) are holomorphic in some neighborhood of and their values 
are the same in the sequence of points having as an accumulation 
point, we obtain W\(t) = W2(t). Thus any solution of Eq. (fT9|) mero- 
morphic in the entire complex plane and holomorphic in a neighbor- 
hood of zero has the form 1)22(1 . 

This is "nearly" Statement 2' from Section [21 "Nearly" because 
the initial point to cannot be arbitrary, but should satisfy the condi- 
tion |io| < R and this R depends on the solution. However it per- 
mits us to find a superposition formula. Indeed, let Wi(t), . . . , W m (t) 
be particular solutions meromorphic in the entire complex plane and 
holomorphic in a neighborhood of zero. For each Wj we can find its 
own constant Ri and take R = min(i?i, . . . , R m ). 

Let us choose an initial parameter value to such that |to| < R- 
It follows that Wx(t), . . . , W m {t) have the form Wi(t) = g(t) ■ Wi(t ), 
where g(t) is the same for all i. As in Section [21 we have the system 
of equations for g(t). If m is sufficiently large and Wj(i) satisfy some 
additional independence conditions necessary for g(t) to be expressed 
in terms of Wi(t) and Wi(to), we have 

g(t) = F(ui(t), . . . ,u m (t);ui(t ), . . .,u m (t )). 

It follows that the general solution holomorphic in the disc \t\ < \to\+e 
has the form 

u(t) = g(t) -u = F(ui(t), . . . ,u m (t);ui(t ), . . . ,u m (t )) -u(t ), (23) 

the initial condition u(to) plays a role of the arbitrary constant in the 
superposition formula (|23|) . Since the function F does not depend on 
the choice of to, we can consider a general solution holomorphic in a 
disc of radius less then R. If we decrease R and repeat our construction 
of the superposition formula, the resulting superposition formula will 
be the same. This finishes the proof. □ 

As we already mentioned in Section [21 the resulting superposition 
formula is the same as in the case of usual Riccati equation © . These 
superposition formulas can be found in 0IE1- 
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4.2 Superposition formulas in the case q = 1, 



By rescaling t = ht' we can transform Eq. (|19|) into the same equation 
with q = 1, h = 1. 

This case is very different from the case |^| 7^ 1. It is difficult to find 
a natural class of solutions defined in the entire complex plane. Since 
the set {to + nh, n € N} has no accumulation point, the uniqueness 
theorem for holomorphic functions cannot be applied. In this situation 
we have, for example, a uniqueness theorem for entire functions of 
growth order 1 and of normal type, but it is very unnatural. We 
already mentioned in the beginning of Subsection l4.1l that solutions of 
the Riccati equations are in general meromorphic. But the situation 
becomes very natural if we consider functions defined only on Z, i. e. 
in integer points. 

Theorem 4 The discrete Riccati equation 



W{n+l)-W{n) = A{n)+B{n)W{n)+W(n+l)C{n)+W(n+l)D{n)W{n) 



are non- degenerate for all n has a superposition formula in the class 
of solutions defined on all Z. 

Proof. Let us consider the formula (|2U|) as a difference equation for 
the N x iV-matrix g(t). We can rewrite it as 



Let us also fix some uq and consider the initial problem g(no) = I. 
Since the matrix (|25|) is non-degenerate, 



(24) 



such that the matrices 




(25) 





is uniquely defined for all n and g(n) G GL(N). 
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Let us consider a solution W\{ri) defined on all Z. Let us use the 
same np and let Wo be the initial condition Wp = Wi(no). We can 
consider a function 



W 2 {n) = g{n) ■ Wp = (M(n)W + N(n))(P(n)Wp + Q(n))" 1 . (26) 



It is also a solution of Eq. (|24|) with the initial condition W^no) = Wo- 
Eq. (|24|1 can be rewritten in the following forms: 



W(n + 1) [I - C(n) - D{n)W{n)] = A{n) + B{n)W{n) + W(n), 

\I + B{n) + W{n + l)D(n)]W{n) = W(n + 1) - A(n) - W{n + l)C(n). 

It follows that if for all n we have 

det[/-C(n)-D(n)Wi(n)] / 0, det[J+B(n-l)+Wi(n)D(n-l)] / 0, 



then for all n we can find T4 /r i(n) and Wviri) starting from the initial 
condition Wo = W\{no) = ^(no) and we can see that Wi{n) = 
W2 (n) . We see that all the solutions satisfying the conditions (|27|) are 
of the form (|26|) , This implies in the same way as in the proof of 
Theorem |31 that equation (|24|) has a superposition formula in the class 
of solutions defined on all Z and satisfying (|27|) , Moreover, a general 
solution satisfies conditions (|27|) since they represent only a countable 
number of inequalities (labeled by n). □. 

In the case of the differential Riccati equations we have a super- 
position formula in the class of C 1 (IR)-solutions. However, for some 
particular choice of the constants this formula gives us solutions out- 
side this class, for example with poles. The same situation arises in the 
case of equation (|24fl . We have a superposition formula in the class 
of solutions defined on all Z, but for some particular choices of the 
constants this formula gives us solutions which are singular in some 
points. 

For example, let us consider 



(27) 




This curve gives us solutions 



w(t) 



(u Q - l)t 2 + Up 
(up - l)t 2 + 1 
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of the equation 

w(t) = -2t + Atw(t) - 2tw(t) 2 

with initial condition w(0) = uq. For uq > 1 these solutions belong 
to C (R) and we can write a superposition formula. If we take three 
particular solutions Wq,w%,W2 corresponding to Uq = l,Uo = 2 and 
uq = 3, respectively, we obtain from (J2J) the formula for the general 
solution belonging to C 1 (IR) : 

2t 2 + 4-c 

Wit) = K . 

w 2t 2 + 2 - c 
If, for example, c = 4, we obtain a solution 



w(t) 



t 2 



t 2 -r 



which is singular at t = ±1. 

In the difference case the same curve gives us solutions 

(n - l)n 2 + uq 

win) = — r — k 

V ' (u -l)n 2 + l 

of the equation 

w(n+l)-w(n) = -2n-l+(2n+l)w(n)+(2n+l)w(n+l)-(2n+l)w(n+l)w(n). 
We obtain the same superposition formula 

2n 2 + 4 - c 



w(n) 



2n 2 + 2-c 



valid for the solutions defined on all Z. It also gives us solutions with 
singularities. For example for c = 4 we have 



n 2 



win) 



n 2 



which is singular at n = dbl. We see that the superposition formula can 
be used to avoid singularities in numerical computations, as already 
pointed out in Ref. [T3] . 
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5 Conclusions 



We have constructed very natural discrete matrix Riccati equations 
with superposition formulas. This was done for a family of difference 
operators general enough to include both (/-derivatives and standard 
discrete derivatives. 

The key idea was considering superposition formulas in an appro- 
priate class of solutions. Which class depends on the type of difference 
operator used. It turns out that in the case of the g-derivative it is 
natural to consider solutions meromorphic in the entire complex plane 
and holomorphic in a neighborhood of zero, and in the case of the stan- 
dard discrete derivative, it is natural to consider solutions defined only 
in integer points. 

The discretization of Riccati equations of this article should be 
compared with a previously proposed discretization of linearizable 
equations, including the Riccati one ^1 E3- There the curve in a 
Lie algebra, i. e. of Eq. Q, was replaced by a curve in the corre- 
sponding Lie group G and this gave a linearizable mapping. E. g. for 
the Riccati equation the corresponding discretization would be 



w(t + h) = [g u {t)w(t) + g 12 {t)}[g 21 {t)w{t) + g 22 (t)}-\ (28) 



Here h is, on one hand, a lattice spacing, on the other hand a group 
parameter. Thus, the whole class of Riccati equations is replaced 
by the above class of homographic mappings. If equation ()28|) has 
constant coefficients, then the relation between the Riccati mapping 
and the Riccati equation is simple: 



i. e.: g u = 1 + hb + . . . , g u = ha + . . . , g 2 i = -hd + . . . , g 22 = 
1 — he + . . . , and for h — > (|2*H|) reduces to 



For variable coefficients the relation between a specific Riccati map- 
ping and a specific Riccati equation is difficult to establish (i. e. it is 
necessary to integrate the Riccati equation explicitly). For all details 
we refer to the original articles |14| 115] . We also mention that discrete 
integrable mappings of the type have been used to discretize the 
Pinney equation ]16j . 




w(t) = a + (b + c)w(t) + dw 2 (t). 
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The relation between the discrete and continuous equations of 
this article, on the other hand, is quite straightforward. Given the 
ODE (j3J), we give the discrete equation (fT§]) explicitly, involving the 
same functions A(t), . . . , D(t). Thus equation (|T§|) can serve as the 
basis for numerical approximations for differential equations. This 
would complement a previous use of superposition formulas in numer- 
ical analysis |13| . 

We would like to comment on the relation between the discrete 
Riccati equations and the Runge-Kutta method. Statement 2 of Sec- 
tionElstates that in the class of C 1 (IR)-solutions the differential Riccati 
equation (j5J) is equivalent to the linear matrix equation 

««> = ( ^1W> <») 



■D(t) -C(t) 

with the initial condition g(to) = I. Applying the Runge-Kutta 
method to Eq. (|29|). we obtain the difference equation 

= -^)) 9(t) 

which is in a similar way equivalent to the difference Riccati equation 
defined by the operator U\ t h- As it was explained before this will also 
give us singular solutions. 

Equivalently we can say that all C 1 (M)-solutions of the Riccati 
equation © have the form W(t) = X(t)Y(t)- 1 , where X(t) and Y(t) 
are solutions of the linear system 

d_ ( X{t) \ _ ( B(t) A(t) \ ( X(t) \ 

dt \ Y(t) ) \ -D(t) -C(t) ) \ Y(t) ) ■ [6U) 

Applying the Runge-Kutta method to Eq. (|3fl|) gives us the equation 

/ X(t) \ _ ( B(t) A(t) \ ( X{t) \ 
Ul ' h { Y(t) ) ~ { -D(t) -C(t) ) { Y(t) ) ■ ^ 

Introducing the inhomogeneous coordinates (the matrix elements of 
W = XY' 1 ) we transform (|31|) into the discrete Riccati equation ()19|) 
with q = 1, h ^ 0. 

In other words, the /i-discretization of the matrix Riccati equa- 
tion, presented in this article, not only preserves the superposition 
formula. It is also equivalent to applying the Runge-Kutta method to 
the associated linear system (|30|) . 
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